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ABSTRACT 



We present new near-IR H2, CO J=2-l, and CO J = 3-2 observations to study 
outflows in the massive star forming region IRAS 05358+3543. The Canada-France- 
Hawaii Telescope H2 images and James Clerk Maxwell Telescope CO data cubes of the 
IRAS 05358 region reveal several new outflows, most of which emerge from the dense 
cluster of sub-mm cores associated with the Sh 2-233IR NE cluster to the northeast of 
IRAS 05358. We used Apache Point Observatory (APO) JHK spectra to determine line 
of sight velocities of the outflowing material. Analysis of archival VLA cm continuum 
data and previously published VLBI observations reveal a massive star binary as a 
probable source of one or two of the outflows. We have identified probable sources for 
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6 outflows and candidate counterffows for 7 out of a total of 11 seen to be originating 
from the IRAS 05358 clusters. We classify the clumps within Sh 2-233IR NE as an 
early protocluster and Sh 2-233IR SW as a young cluster, and conclude that the outflow 
energy injection rate approximately matches the turbulent decay rate in Sh 2-233IR NE. 

Subject headings: ISM: jets and outflows — ISM: Herbig-Haro objects — ISM: kine- 
matics and dynamics — ISM: individual: S233IR — stars: formation 



1. Introduction 



Collimated, bipolar outflows accompany the birth of young stars from the earliest stages of 
star formation to the end of their accretion phase (e.g. iReipurth &: Ballvll200ll ). While the birth 
of isolated low-mass stars is becoming well understood, the formation of massive stars (> WM & ) 
and clusters remains a topic of intense study. Obse rvations show that moderate to high- mass stars 
tend to form in dense clusters (jLada fc Ladall2003l ). In a clustered environment, the dynamics of 
the gas and stars can profoundly impact both accretion and mass- loss processes. Feedback from 
these massive clusters may play a significant role in momentum injection and turbulence driving in 
the interstellar medium. 

Outflows from massive stars are less studied than those from low mass stars largely because 
massive stars accrete most of their mass while deeply embedded. Therefore, unlike low mass young 
stars that are accessible in the optical, massive stellar outflows can only be seen at infrared and 
longer wavelengths. Direct evidence for je ts from massive young stellar objects (YSOs) from H2 



or op tical emission is generally lacking (e.g. lAlvarez Hoard 120051 : iKumar et al.l 120021 : 1 Wang et al 



2003), although there is evidence that massive stars are the sources of collimated molecular out- 
flows from millimeter observations (e.g. iBeuther et al.ll2002bl ). Outflows from massive stars may 
allow accretion to contin ue after their radiation pressure would otherwise halt accretion in a spher- 
ically symmetric system (jKrumholz et al.l 12009 ) . They therefore represent a crucial component in 
understanding how stars above ~10 Mq can form. 

IRAS 05358 is a double cluster of embedded infrared sources located at a distance of 1.8 kpc 
in the Auriga molecular cloud complex (jHever et al.lll9 96) associated with the HII regions Sh-2 231 
through 235 at Galactic coordinates around I, b = 173.48,-1-2.45 in the Perseus arm. Sh 2-233IR NE 
is the collection of highly obscured and mm-bright sources slightly northeast of Sh 2-233IR SW, 
which is the location of the IRAS 05358+3543 point source and the optical emission nebula (see 
Figure [T]). The IRAS source is probably a blend of the three brightest infrared objects in the MSX 
A-band and MIPS 24 [im images, which are located at Sh 2-233IR NE, IR 41, and IR 6. For the 
purpose of this paper,the whole complex including both sources is referred toas IRAS 05358, and 
otherwise refer to individual objects specifically. 



Early observations revealed the presence of OH (jWouterloot et al 



1993), H 2 (jScalise et al 
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1989 : iHenning et al.lll992l ). and methanol (|Mentenlll99ll ) masers about an arcminute northeast of 
the IRAS source, indicating that massive stars are likely pr esent at that locat i on. Near infrared 



observations revealed the presence of two embedded clusters (jPorras et al.ll2000l ; iJiang et al.ll2001 
labeled Sh 2-233IR SW for the southwestern cluster associated with the IRAS source, and Sh 2- 
233IR NE f or the northeastern cluster located near the OH, H2O, and CH3OH masers. Stars 
identified in IPorras et al.1 (200C) are ref erred to by the desig nation "IR (number)" corresponding 
to the catalog number in that paper. IPorras et al.l (J200CJ) also included scanning Fabry-Perot 
velocity measurement s of the inner ~ 1 ' . CO observations revea l ed br oa d line wings i n dicat ive of 
a molecular outflow (ICasoli et al.l Il986l ; IShepherd Churchwelll Il996l ) . iKumar et al.l (|2002l ) and 
Khanzadvan et al.l (|2004l ) presented narrow band images of 2.12 [im H2 emission that reveled the 
presence of multiple outflows. Interferometric imaging of CO and SiO confirmed t he presence of 



at lea st three flows emerging from the north east cluster ce n tered o n the masers (jBeuther et al 
2002al ) having a total mass of about 20 M . beuther et all (1200231) also presen ted MAMBO 1.2 



mm maps and a mass estimate of 610 M & for the whole region. IWilliams et al.l (|2004l ) presented 
SCUBA maps and mass estimates of th e clusters of 195/126M for Sh 2-233IR NE and 24/12 M Q 
for Sh 2-233IR SW (850 ^m/450 fim). IZinchenko et al.1 (119971 ) measured the dense gas properties 



using the NH3 (1,1) and (2,2) lines. They measure a mean density n ~ io 3 - 60 cm , temperature 
26. 5K, and a mass of 600 Mq. The total lumino sity of the two clust ers is about 6300 Lq, indicating 
that the region is giving birth to massive stars (jPorras et al.l l2000). 



Millimeter wavelength interferometry with arcsecond angular resolution has reve aled a compact 



2002a. 


2007; 


Leurini et al. 


2007) 



Source mmla is associated with a cm continuum point source and will be discussed in detail below. 

IRAS 05358 has previously been o bserved at low spati al resolution in the J=2-l and J=3-2 tran- 
sitions with the Kosma 3m telescope (|Mao fc Zend |2004| ) . While the general prese nce of outflows 



was rec ognized and a total mass estimated, the specific outflows were not resolved. iBeuther et al 



(|2002al ) observe d the CO J=6-5, J=2-l, and J=l-0 transitions at moderate resolution in the inner 
few arcminutes. iThomas &: Fullerl ( 20081 ) observed C 17 in the J=2-l and J=3-2 transitions with 
a single pointing using the JCMT. 



2. Observations 

A collection of data acquired by the authors and from publicly available archives is presented. 
An overview of the data is presented in figure [TJ The goal was to develop a complete picture of 
the outflows in IRAS 05358 and their probable sources. CO data were acquired to estimate the 
total outflowing mass and to identify outflowing molecular material unassociated with H2 shocks. 
Archival Spitzer IRAC and MIPS 24 [im data were used to identify probable YSOs as candidate 
outflow sources. Near-infrared spectra were acquired primarily to determine H2 kinematics and 
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develop a 3D picture of the region. Optical spectra were acquired to attempt to identify stellar 
types in the unobscured Sh 2-233IR SW region. Finally, archival VLA data were used to acquire 
better constraints on the position and physical properties of the known ultracompact HII (UCHII) 
region, and to detect or set limits on other UCHIIs. 

2.1. Sub-millimeter Observations 

The 345 GHz J = 3-2 rotational transition of CO was observed with the James Clerk Maxwell 
Telescope (JCMT) on 4 January, 2008 with the 16 element (14 functional) HARP-B heterodyne 
focal plane array. Two 12' x 10' raster scans in R.A. and Dec. were taken with orthogonal 
orientations to assure complete coverage in the region of interest; this resulted in a useable field 
11.7' x 11.3' with higher noise along the edges. The beam size at 345 GHz is about 15". 

Observations were conducted during grade 3 conditions with the 225 GHz zenith optical depth 
of the atmosphere r ~ 0.1. A channel width of 488 kHz corresponding to 0.423 km s _1 was used. 
The maps required a total of 1 hour to acquire and resulted in an effective integration time of 4.6 
seconds per pixel (there are 12,000 6 x 6" pixels in the final grid), resulting in a noise per pixel of 
0.36 K km s -1 . 

The optical depth and telescope efficiency corrections were applied by the JCMT pipeline to 
convert the recorded antenna temperatures to the corrected antenna 0. An additional main-beam 
correction has been applied, 

Ta* 
'j - ■ i 

Vmb 

where rj m b was measured by observing Mars to be ~ 0.60 at 345 GHz. Emission in the sidelobes is 
expected to be small at the outflow velocities. 

On September 25 and November 15, 2008 the CO, 13 CO, and C 18 J=2-l transitions were ob- 
served in the central 3' of IRAS 05358. The beamsize at 220 GHz is about 23". The sideband config- 
uration used also includes the 
SO 56 — 45 and 13 CS 5-4 transitions. Conditions during these observations were grade 5 (r ~ 
0.24 — 0.28) and therefore too poor to use the HARP instrument, but acceptable for the A3 detec- 
tor. 

Data reduction used the Starlink package following the standard routines recommended by the 
JCMT support scientists El- The CO 3-2 data cube was extracted over a velocity range from -50 to 
10 km s _1 LSR and spectral baselines were fit over the velocity range -50 to -40 and to 10 km s _1 
and subtracted. The data were re-gridded into 6" pixels and 2 pixel Gaussian smoothing was used 



'See 

http://docs.jach.hawaii.edu/JCMT/OVERVIEW/tel_overview/ for a discussion of JCMT parameters 
2 http: //www. jach.hawaii . edu/JCMT/spectral_line/data_reduction/acsisdr/ 
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to fill in the gaps left by the two bad detectors in the 4x4 array. The data cube was cropped to 
remove undersampled edges which have high noise and bad baselines. The beam efficiency was 0.68 
at 230 GHz. 

The A3 data cubes were extracted over the velocity range -60 to 20 km s _1 and baselines were 
calculated over -60 to -40 and to 20 km s" 1 . The data was gridded into 10" pixels with 2 pixel 
gaussian smoothing to reduce sub-resolution noise variations. 



2.2. Spitzer 

Spitzer IRA C bands 1 to 4 and MIPS band 1 data were retrieved from the Spitzer Science 



Center archive. iQiu et al.l (|2008l ) acquired the data as part of a study of many high-mass star 



forming regions; they identified YSO candidates based on IRAC colors. T he version 18 po st-BCD 
data products were used to produce images and photometric catalog from lQiu et al.l (|2008l ). which 
was made from a more carefully-reduced data set, was used for SED analysis. 



2.3. Near-IR images 

Near-infrared data were acquired using the Wide-field Infrared Camera (WIRCam) on the 
Canada- France-Hawaii Telescope (CFHT) on Mauna Kea. The field of view is 20'x20' and pixel 
scale 0.3". Data were acquired on November 18, 19 and December 20, 2005. The seeing was 0.5-0.7" 
during the observations. A 0.032 \xm wide filter centered at 2.122 jim was used to take images of 
the H2 S(l) 1-0 rovibrational transition. Each H2 exposure was 58 seconds, and dithered images 
were taken for a total exposure time of 1755 seconds. The data were reduced with the WIRCam 
pipeline. 



2.4. Near-IR spectra 

Near-infrared spectra were acquired using the TripleSpec instrument at Apache Point Obser- 
vatory. TripleSpec simultaneously acquires J, H, and K band spectra over a 42" long slit. A slit 
width of l.l" with an approximate spectral resolution A/AA = 2700 was used. 

Observations were taken on the nights of December 2, 2008 and January 7, 2009. Data on 
December 2 were taken in an ABBA nod pattern, but because of the need to observe extended 
structure across the slit a stare strategy was selected on January 7. 

The data were reduced using the twodspec package in IRAF. HD31135, an AO star, was used 
as a flux calibrator. Wavelength calibration was performed using night sky lines. Lines filling the 
slit were subtracted to remove atmospheric emission lines. Telluric absorption correction was not 
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performed, but telluric absorption is considered in the analysis. 

The transformations from the observed geocentric reference frame to vlsr were computed to 
be 0.78 km s -1 on Dec 2 and 19.74 km s -1 on Jan 8. 



2.5. Optical Spectra 

Optical spectra were acquired using the Double Imaging Spectrograph instrument at APO. 
The high-resolution red and blue gratings were centered at 6564 A and 5007 A with a coverage of 
about 1200 angstroms and resolution A/AA rs 5000. Sets of three 900s exposures and three 200s 
exposures were acquired on the targets and on the spectrophotometric calibrator G191-b2b with a 
1.5" slit. Observations were taken on the night of January 17, 2009 under clear conditions. 

Optical spectra were also reduced using the twodspec package in IRAF. Wavelength cali- 
bration was done with HeNeAr lamps and night sky lines in the red band, and HeNeAr lamps in 
the blue band. Lines filling the slit were subtracted to remove atmospheric lines, though some 
astrophysical lines also filled the slit and these were measured before background subtraction. The 
vlsr correction for this date was 24.4 km s -1 . 



2.6. Optical imaging 

CCD images images were obtained on the nights of 14 and 15 September 2009 NOAO Mosaic 
1 Camera at the f/3.1 prime focus of the 4 meter Mayal telescope atthe Kitt Peak National Obser- 
vatory (KPNO). The Mosaic 1 camera is a 8192x8192 pixel array (consisting of eight 2048x4096 
pixel CCD chips) with a pixel scale of 0.26" pixel -1 and a field of view 35.4' on a side. Narrow-band 
filters centered on 6569A and 6730A both with a FWHM of 80A were use to obtain Ha and [SII] 
images. An SDSS i' filter which is centered on 7732A with a FWHM of 1548Awas used for contin- 
uum imaging. A set of five dithered 600 second exposures were obtained in Ha and [SII] using the 
standard MOSDITHER pattern to eliminate cosmic rays and the gaps between the individual chips 
in Mosaic. A dithered set of five 180 second exposures were obtained in the in the broad-band SDSS 
i-band filter to discriminate between Ha, [SII], and continuum emission. Images were reduced in 



the standard manner by the NOAO Mosaic reduction pipeline (jValdes fc S waters! 120071 ). 



2.7. VLA data 



VLA archival data from projects AR482, AR513, AS831, and AM697 were re-reduced to 
perform a deeper search for UCHII regions and aquire more data po i nts on the known UCHII's 
SED. Data from AR482 were previously published in iBeuther et al.l (|2007l ). the other data are 
unpublished. The data were reduced using the VLA pipeline in AIPS (vlarun). The observations 
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Fig. 1.— The CFHT H 2 (bold), CO 3-2 HARP (thin), and CO 2-1 A3 (dashed) fields overlaid on 
the KPNO Ha mosaic with selected objects identified by their SIMBAD names. Sh 2-233IR SW is 
coincident with IRAS 05358+3543. 
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used and sensitivities and beam sizes achieved are listed in Tables [T] and There appeared to be 
calibration errors in the AR482 observations (the phase calibrator was 2-3 times brighter than in 
all other observations) and this data were therefore not used in the final analysis, but it produced 
consistent pointing results. 



3. Results 



3.1. Near Infrared Imaging: Outflows and Stars 



Eleven distinct outflows have been identified in IRAS 05358 in the images. Outflows are 
identified from a combination of J=3- 2 CO data, shock excited H2 emission, and published inter- 
ferometric maps (jBeuther et al.ll2002al ). Suspected CO outflows were identified by the presence of 
wings on the CO J=3-2 emission lines that extended beyond the typical velocity range of emission 
associated with the li ne core. The single dish data were compared to the interferometric maps of 
Beuther et al.l (|2002al ). The CFHT H2 image was then used to search for shock-excited emission 



associated with the outflow lobes. 

Figure [5] shows the H2 S(l) 1-0 2.1218 fim (a rovibrational transition in the electronic ground 
state from the v = 1, J = 3 to the v = 0, J = 1 state) emission in t he vicinity of IRAS 5358 
with outflows and possible outflow sources labeled. The mm cores from lBeuther et al.l (|2002al ) are 
identified by red squares. 

The flow vectors in figure [2] were chosen on the basis of the H2 bow shock morphologies 
and orientations of c hains of H2 features, association with arcsecond-scale CO features on the 
Beuther et al.l (|2002al ) Figure 8 CO map, and/or association with lobes of Doppler-shifted CO 
emission in the CO 3-2 data. The color of the vector indicates the suspected Doppler shift; red and 



Table 1. VLA Observation Program Names, Dates, and Times 



VLA Observation Time Array Band Fluxcal Phase cal Phase cal 

Observation Date on Percent 

Name Source Uncertainty 



AR482 


August 2 2001 


2580s 


B 


X 


3c286 


0555+398 


22 


AR513 


June 21 2003 


7770s 


A 


X 


3c286 


0555+398 


0.8 


AS831 


February 26 2005 


2640s 


B 


X 


3c286 


0555+398 


0.7 


AS831 


August 5 2005 


2660s 


C 


X 


3c286 


0555+398 


0.3 


AS831 


May 11 2006 


2610s 


A 


X 


3c286 


0555+398 


3.0 


AL704 


August 7 2007 


6423s 


A 


Q 


3c273 


0555+398 


18 


AL704 


September 1 2007 


6423s 


A 


Q 


3c273 


0555+398 


13 


AM697 


November 26 2001 


2880s 


D 


Q 


3c286 


0555+398 


2.2 


AM697 


November 28 2001 


1530s 


D 


K 


3c286 


0555+398 


2.1 


AM697 


November 28 2001 


1530s 


D 


u 


3c286 


0555+398 


5.8 
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Fig. 2. — The outflows described in section [37T1 overlaid on the CFHT H2 image. Numbers followed 
by r and b (red and blue), n and s (north and south), or e and w (east and west) are thought to 
be counterflows. Red and blue vectors indicate red and blue doppler shifts. Green vectors indicate 
where the doppler shift i s ambiguous o r canno t be determined. Magenta circles are Spitzer 24^m 
sources. Red squares are lBeuther et al.l (|2002al ) mm sources (from left to right, mml, mm2, mm3). 
The blue diamond is a YSO candidate detected only in IRAC bands. The length o f the vectors 



corre sponds to the approximate length of the outflows. Source 1 and 6 correspond to lPorras et al 



(|2000l ) IR 6 and IR 41 respectively, and they are discussed under these names in sections 13.11 The 
bows of Outflow In and 4n are detected in Ha and [S II] emission and are therefore as identified 
as Herbig-Haro objects HH 993 and 994 respectively. 
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Fig. 3. — CO contours integrated from vlsr = -13 to -4 km s 1 (red) and -34 to -21 km s 1 (blue) 
at levels of 2,4,6,10,20,30,40,50 K km s _1 overlaid on the H 2 image. Specific outflows are labeled 
in Figure [2] on the same scale. 



blue correspond to red and blueshifts and green vectors indicate that the Doppler shift is uncertain. 



IRAS 05358 outflow 1: The most pro minent flow in H2 is associated with the bright bow- 
shocks Nl and N6 (jKhanzadyan et alj|2004l) locat ed towards P A 345° and 170° res pectively from 
the sub-mm source mmslb (jBeuther et all l2002al ) ■ This flow, iBeuther et al.1 (|2002al ) outflow A, is 
associated with redshifted and blueshifted CO emission. The northern shock is seen in Ha and [S 
II] emission (figure Hb) and is given a Herbig-Haro designation HH 993. 

This flow is indicated by oppositely dir e cted green vectors f rom the vicinity of smml, 2, and 



3. It is listed as "Jet 1" in lQiu et al.1 (120081 ) 



behind the bow shock as a Mach disk. In the Beuther et al 



Kumar et al 



(120021) identified the knot immediately 
(|2002al ) interferometric maps, the north 



flow contains redshifted features and the south flow contains primarily blueshifted features. There 
are also blueshifted CO features to the west of the H2 knots that are probably part of a different 
flow that is not seen in H2 emission. 

The velocity of the flow as measured from H2 emssion is blueshifted as much as 80 km s _1 (LSR), 
but one component is blueshifted only 14 km s _1 (see table [3]), which is consistent with the cloud 
velocity. A redshifted SiO lobe is present in the south counterflow. The presence of Ha, [S II], and 
[O III] emission in the north shock and corresponding nondetections in the south shock suggest that 
there is substantially greater extinction towards the south knot. While the velocities in three of the 
four apertures picked along the TripleSpec slit are blueshifted, there are also knots with velocities 
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Fig. 4. — (a) H2 image with SO 56 — 4s pea k flux contours at 0.5-1.4 K in intervals of 0.15 K 
overlaid. With a critical density ~ 3.5 x 10 6 ( Schoier et al.ll2005l ). this transition is a dense gas 
tracer, (b) The [S II] image with outflow vectors overlaid. Diffuse emission can be seen at the north 
ends of Outflows 1, 4, and 6 and around the reflection nebula near source IR 41. 



consistent with the cloud velocity. iPorras et al.1 (120001 ) measure the velocity of the counterflow to 
be -17.3 km s , which is consistent with the cloud velocity. Outflow 1 is propagating very nearly 
in the plane of the sky. 



A line connecting the two bow shocks in Outfl ow 1 goes directly thro ugh iBeuther et al.l (|2007l ) 
source mm2a despite the clear association in the IBeuther et al.l (|2002al ) interferometric CO map 
(their Figure 8) with mmla. The currently available data do not clarify which is the source of the 
outflow: while the bent CO outflow appears to trace Outflow 1 back to mmla, there are additional 
parallel CO outflows towards the confused central region that could originate from either mmla or 
mm2a. 

A Spitzer 4.5 \im and 24 jim source is barely detected in H2 2.5' to the north of Outflow 1. It 
is only apparent when the H2 image is smoothed and would have been dismissed as noise except 
for the association with a probably 4.5 /xm extended source. It is labeled 2A[im source 7 in figure 
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[2j It appears to be slightly resolved at 4.5//m, and is therefore likely shocked emission. The object 
may be a protostellar source with an associated outflow, but its proximity to the projected path of 
Outflow 1 suggests that it may be an older outflow knot. 

IRAS 05358 Outflow 2: The second brightest H2 features trace a bipolar flow emerging from 
the immediate vicinity of the s ub-mm cluster a t PA ~ 135° (red lobe) and 315° (blue lobe). It is 
listed as "Jet 2" in Figure 6 of iQiu et al.l (|2008l ). The counterflow probably overlaps in the line of 
sight with the counterflow from Outflow 3. It is shorter on the counterflow side either because it 
has already penetrated the cloud and is no longer impacting any ambient gas or, more likely, it has 
slowly drilled its way out of the molecular cloud and has not been able to propagate as quickly as 
the northwest flow. The H2 velocities measured for these knots are ~ 30 km s _1 blueshifted, or 
marginally blue of the cloud LSR velocity. 



The disk identified in iMinier et al.l (|2000l ) is approximately perpendicular to the measured 
angle of Outflow 2 assuming that mmla is the source of this flow. It is therefore an excellent 
candidate for the outflow source. A diagram of the mmla region is shown in figure [TBI See Section 
for detailed discussion. 



outflow B at PA 



200 



Khanzac 


van et al. 


Jet 3" in 


Oiu et al. 



The two chains of H2 emission indicate that outflows 2 and 3 are distinct. There also 
appears to be a counterflow at a shorter distance from the mm cores similar to counterflow 2. 

Outflows 2 and 3 may be associated with either redshifted or blueshifted features in the 
Beuther et all (|2002a | X CO and SiO maps. High velocity flows with both parities are present near 
both the northwest (jBeuther et al.1 (|2002al ) outflow C) and southeast flow for these jets, but the 
resol ution of the millime ter observations is inadequate to determine which flow is in which direc- 

1 for their knot 4A, which corresponds to 



tion. 



Porras et al.l (|2000l ) measures vlsr = —7.5 km s" 
the blended southeast counterflow of outflows 2 and 3. Their Figure 7 shows a wide line that is 
probably better represented by two or three blended lines, one consistent with the cloud velocity 
and the other(s) redshifted. Since Outflow 2 has a measured blueshift and outflow 3 is signifi- 
cantly fainter, the redshifted counterflow emission is probably associated with Outflow 2 and the 
blueshifted with outflow 3. 

IRAS 05358 outflow 4: The JCMT CO data and H2 images reveal a large outflow lobe consist- 
ing of blue lobes 1 and 4 that form a tongue of blueshifted emission propagating to the northeast at 
PA « 20° (Figure [2]) from the cluster of sub-mm cores. A faint chain of H2 features runs along the 
axis of the CO tongue and terminates in a bright H2 bow shock located at the northern edge of [2j 
Several H2 knots lie along the expected counterflow direction, but that portion of the field contains 
multiple outflows and is highly confused. If the counterflow is symmetric with the northeast knot, 
it extends 2.1 parsecs on the sky. 



The bow shock of Outflow 4 is seen in the HII and [S II] images, implying that the extinction 
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is much lower than in the cluster. Two apertures placed along the bow shock reveal that it is 
blueshifted about 70km s _1 and may be extincted by as little as Ay ~ .5. It is designated HH 994. 

IRAS 05358 outflow 5: Figure[2]shows a bright chain of H2 knots and bo w shocks starting abou t 
10" west of mm3 and propagating south at PA ~ 190°. The SiO maps of iBeuther et all (|2002al ) 
show a tongue of blueshifted emission along this chain (their Outflow C). The outflow projects 
back to H 13 CO + source 3, which is also a weak mm source. A lack of obvious counterflow and the 
possibility that the knots identified with Outflow 5 could be associated with a number of different 
crossing flows makes this identification very tentative. Higher spatial resolution observations will 
be required to determine the association of this outflow. 

IRAS 05358 outflow 6: The fourth brightest source in the Spitzer 24/im data is loc a ted a t 
J(2000) = 05:39:08.5, +35:46:38 (source 5 in the IRAS 05358 section of the loiu et al.1 (hoosh 
catalog, referred to in table [3] as Q5) in the middle of the molecular ridge that extends from IRAS 
05358 towards the northwest (24^m object 4 in figure El). The sta r is located at the northwest end of 
the tongue of 1.2mm emission mapped by (jBeuther et al.ll2002al ) with the MAMBO instrument on 
the IRAM telescope. This part of the cloud is also seen in silhouette against brighter surrounding 
emission at 8/im. At wavelengths below 2fim, it is fainte r than 14-th magnitude and therefore is 
not listed in the 2MASS catalog, and it is not detected in lYanl (|2009l ) down to 19th magnitude in 
K. 

Spitzer data indicates very red colors between 3.6 and 70 /im, indicating th at this object is 



likely to be a Class I protostar. The SED is fit using the online tool provided by iRobitaille et al 



(|2007l ). Unfortunately, a wide variety of parameters all achieved equally good fits, so no conclusions 
are drawn about the stellar mass or other very uncertain parameters. However, the top models all 
had Ay > 20 and many in the range 30-50, indicating that the line of sight is probably through a 
thick envelope or disk towards this source. 

This source lies at the base of the tongue of blueshifted CO 3-2 emissi on that extends northwes t 
of IRAS 05358 at PA « 345° and has mass ~ .5M . A pair of H 2 features, Khanzadyan et al.1 (|2004l ) 
N12A and N12B are located 30 and 55" from the susp ected YSO, forming a chain along the axis of 
the blueshifted CO tongue. Khanzadyan et al. ( 2004! ) H2 knot N3F lies along the flow axis in the 
redshifted direction. 



IRAS 05358 outflow 7: The 20" long chain of H 2 knots labeled iKhanzadvan et al.1 (120041 ) Nil 
appears to trace part of a jet at PA w 345° tha t propagated parallel to outflow 6 about 20" to the 
east. The northwest portion of Outflow C in the lBeuther et al.1 (|2002al ) SiO map is in approximately 
the same direction as Outflow 7, and it may represent a redshifted counterflow to the northwest- 
pointing H 2 knots. The jet axis passes within a few arc-seconds of a faint and red YSO located at 
J(2000) = 05 39 10.0, +35 46 27 (blue diamond in figure [2] about 35" south of the southern end 
of the H2 feature). It may be a 24^m source but is lost in the PSF of the bright so urce at the 
center of Sh 2-233IR NE. This object is also undetected down to 19th magnitude in the lYanl ((2009) 
K-band image. 
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IRAS 05358 outflow 8: A prominent jet-like H2 feature protrudes from the vicinity of Sh 2- 
233IR SW at PA « 335° and ends in bright knot N9. The feature N5B is is located just outside the 
ring of H2 emission that surrounds the IRAS source at the base of the jet. Towards the southeast, 
knot N6 is located opposite knot N9 with respect to the southwest cluster. IR 41, the Ha emission 
source, labeled 24/^m source 6 in figure [21 is probably the source of this outflow. 

IRAS 05358 outflow 9: In the Spitzer and K s images, an infrared reflection nebula opens 
towards the southwest at PA ~ 245° and points towards a blueshifted CO region. The reflection 
nebula is also seen in Ha. It is likely that the CO emission in CO Region 1 (table [2]) traces a fossil 
cavity whose walls provide the scattering surface of the reflection nebula. 

IRAS 05358 outflow 10 and IR 6: A b r ight H ? filament protrudes at PA w 15° towards the 



northeast of IR 6 (24//m source 1. lOiu et al.1 (|2008l ) source 8). The star is the third brightest 24/im 



source in the IRAS 05358 region. Since it is visible at visual wavelengths, it is not heavily embedded. 
Its Ha emission and association with an outflow lobe and H2 emission suggest that it is a moderate 
mass Herbig AeBe star associated with the IRAS 05358 complex. The optical spectrum confirms 
this hypothesis: the star has Ha absorption wings on either side of a very bright, asymmetric Ha 
emission profile (see section 13 . 5[) . 



IR 6 is seen to be the source of Outflow 10. Data for this source is available from ~0.45- 



24/xm, so the iRobitaille et al.1 (|2007l ) spectral fitter puts strong constraints on the star's mass and 
luminosity. The measured mass and luminosity are M = 4.5 ± 0.5 M and L = 10 2 ' 3± ' 25 -Lq, 
parameters consistent with a B7V (±1 spectral class) main sequence star. The range of ages in the 
models covers 10 4 — 10 7 years but favors stars in the range 10 5 — 10 6 years. 

While there is a small clump of redshifted CO emission to the northeast of the object, the H2 
spectrum shows that the north flow is blueshifted vlsr ~ —40km s _1 , and the lack of a visible 
counterflow suggests that the counterflow may be masked behind an additional extincting medium. 
The counterflow drawn in figure [2] is not seen in emission but is identified as a probable location 
for a counterflow because of the confident association of outflow lOn with source IR 6. 

IRAS 05358 outflow 11: A chain of H2 knots is seen at 2.12^m and in the Spitzer 4.5/im 
image. They trace back to either IR 78 or 24//m source 4. There is a tongue of redshifted CO 3-2 
emission in the same direction as this flow that suggests it may be redshifted. 

IR 41' There is an arc- like H2 emission feature surrounding the Ha emission line star IR 
41. This implies that the star is probably a late B-type star with too little Lyman continuum 
emission to generate a photon-dominated region (PDR) but enough soft UV to excite H2. From 
the measured Ha and nondetection of H/3 at the star's location down to a 5-a l imit of Ixl0~ ] ' 



ergs 1 cm 2 A 1 , a lower limit on the extinction column Ay = 15 is derived. The IRobitaille et al 
(|2007l ) fitter yields a mass estimate of 7.4±O.6M and luminosity L = 1O 2 ' 97± ' 16 L0 among the 222 



best fits out of a grid of 200,000 model SEDs (fits with \ 2 < 5000). The luminosity is very well 
constrained, varying only modestly to L = lO 2 - 99± - 15 L for the 904 best fits (x 2 < 10000), while 
the mass shifts down to 6.5 ± 1.0M Q . The mass estimate may be biased by the lower number of 
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high-mass models computed. The star's mass is most compatible with a main sequence B4V star, 
though its luminosity is closer to a B5V star. The disk mass is constrained to be > 10 3 Mq. The 
age is reasonably well constrained to be T = io 5 - 78±12 for the best 904 models, but is essentially 
unconstrained for the best 222. Similarly, the stellar temperature is entirely unconstrained by the 
fitting process. 

The very high values of x 2 would normally be worrisome, but the x 2 statistic only represents 
statistical error, while the data is dominated by various systematic errors including calibration 
offsets in the optical/NIR and poor resolution in the far-IR. Therefore, it is not possible to find a 
perfect model fit, but still possible to put constraints on the physical properties of the source. 

South of IRAS 05358: There is a symmetric flow with one faint H2 knot and a bright central 
source about 4' south of IRAS 05358. The H 2 knot is at J(2000) = 05:39:15.63 +35:42:13.2. The 
flow has a clear red and blue region as identified in figure [6J the red flow extends from -9 to -14 
km s _1 and the blue from -19 to -23 km s _1 (the outflow is swamped by ambient emission in the 
intermediate velocity range) . The outflow is ~ 2' long, though the probable source identified is not 
directly between the two lobes. The ellipses used are labeled in table [2] as Red S and Blue S. 



3.2. Imaging results: Optical 



Deep [S II] images show that some of the out flows have pierced t hroug h the obscuring dust 
layers or excited extremely bright sulfur emission. iKhanzadyan et al.1 ([2004]) knot Nl at the end 
of Outflow 1 is visible [S II] emission The bow of outflow 4 and the northwest end of outflow 6 are 
detected in [S II]. Only the Outflow 1 and 4 bows are detected in Ha emission, indicating that the 
emission is most likely from shock heating, not external photoionizing radiation. If the shocks were 
externally irradiated, we would expect the emission to be dominated by the recombination lines. 
Because they have been detected in the optical, these two flows can be classified as Herbig-Haro 
objects. 



3.3. CO results 

IRAS 05358 is located at the center of the CO 3-2 integrated velocity maps (Figured]). The 
parent molecular cloud, centered at vlsr = —17.5 km s^ 1 , extends from the southeast towards 
the northwest with the brightest emission coming from the core associated with Sh 2-233IR SW, 
while the highest integrated emission is associated with Sh 2-233IR NE. Sh 2-233IR NE has a 
central velocity of ~ —16.0 km s" 1 from the optically thin C 18 2-1 measurements. Material that 
has been swept up and accelerated by jets and outflows can be seen at velocities vlsr < — 21 
km s _1 and vlsr > ~ 12 km s -1 (Figure [6|). The integrated CO 3-2 map peaks at J(2000) = 
05:39:12.8 +35:45:55, while the highest observed brightness temperature is at J(2000) = 5:39:09.4 
+35:45:12. This offset is discussed in the context of CO isotopologues in section [4 . 4 . 2 1 and in section 
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Fig. 5. — The Ha image with CO contours at redshifted, blueshifted, and middle velocities in 
red, blue, and green respectively. Contours are at 2,4,8,12,20 K km s _1 for the red and blue, and 
20,25,30,40,50,60,70 K km s _1 for the green. Red is integrated from -12 to -4 km s _1 , Blue from 
-31 to -21 km s _1 , and green from -21 to -12 km s _1 . 
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Regions with line wings relative to the ambient cloud within 5' of the northeast cluster were 
assumed to be associated with outflows from the cluster. Further than 5', it is likely that the high 
velocity wings are accelerated by neighboring HII regions (see section [53]) . These line wings were 
integrated over the velocity range -34 to -21 km s _1 (blue) and -12 to 1 km s _1 (red) to acquire 
estimates of the outflowing mass under the assumption that outflowing gas is optically thin. The 
extracted regions are displayed in Figure [6b and measurements in table (2[ The line wings in the 
central arcminute and central 5 arcminutes were measured for comparison with lower resolution 
data and to compute a total outflow mass in the central region. 

The objects in Table[2]labeled CO Region 1, 2, and 3 have uncertain associations with outflows. 
CO Region 1 is tentatively associated with outflow 11. CO region 2 may be associated with Outflow 
3 but is in a highly confused region and may have many contributors. CO region 3 is probably 
associated with outflow 10. In contrast, the associations with outflows 4 and 6/7 are more certain 
because the y are further from th e central region and less confused. Outflow 1 is seen at high 
velocities in 



Beuther et al.1 (|2002al ) interferometer maps. Outflow 9 is selected primarily based on 



CO emission. 



3.4. Near-infrared spectroscopy: Velocities 

The slit positions used and apertures extracted from those slits are displayed in Figure [TUl 
Position-velocity diagrams of the 1-0 S(l) line are displayed in Figure HU Velocity measurements 
are presented in Table [3j 

The near-IR spectrum of Outflow 1 has the largest signal. All of the K-band H2 lines except 
the 2-1 S(0) 2.3556 fim (too weak) and 1-0 S(4) 1.8920 ^m(poor atmospheric transmission) lines 
were detected (see Table [5]). Velocities from gaussian fits to each line are reported. In the central 
portion of Sh 2-233IR NE, outflowing H2 emission at vlsr ~ —30 km s _1 is detected. This material 
may be associated with a line-of-sight flow, or may originate from the base of the already identified 
flows 1-3. In source IR 58, Br7 and He I 2.05835/imare detected, indicating that there is an 
embedded PDR in this source. There is a hint of a second, fainter star adjacent to IR 58. IR 93 
is observed to be a double source in the TripleSpec spectrum, but the spectrum is too weak to do 
any identification. Br7 and possibly He I are detected at fainter levels. 

Table [5] shows the measured line strengths (when detected) of all H2 lines in each aperture. 
The errors listed are statistical errors that do not include the systematics errors introduced by a 
failure to correct for narrow atmospheric absorption lines. 
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Fig. 6. — The JCMT HARP CO J=3-2 map integrated over all velocities with significant emission 
(-34 km s _1 to -4 km s -1 ) shown in gray log scale from to 150 K km s _1 . The elliptical regions 
over which line wings were integrated are shown with blue and red circles corresponding to blue 
and red line wings. The measurements are presented in table [21 
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Fig. 7. — SCUBA 850/um image in linear grayscale from -1 to +10 mJy/beam, with a saturated 
peak of 24 mJy/beam, with 12 CO 2-1 (orange solid, contours at 45,60,85,100,115,130,145 K km s _1 ) 
and 13 CO 2-1 (green dashed, contours at 20,30,40,45 K km s" 1 ) integrated contours. The box shows 
the region plotted in Figure [HI 
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Fig. 8.— CO spectra of Sh 2-233IR NE in 12 CO (blue), 13 CO (green), and C 18 (red). The top-left 
plot is the pixel centered at J(2000) = 5:39:13.67 +35:46:26.0 and each pixel is 10" on a side. The 
region mapped here is shown with a box in Figure UJ Redshifted self-absorption, a possible infall 
tracer, is evident in the 12 CO spectra in the outer pixels. The inner pixels show self-absorption 
only at central velocities: this may be an indication that emission from outflows dominates any 
infall signature, or simply that there is no bulk infall towards Sh 2-233IR NE. 
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Table 2. Measured properties of CO flows 



a Rpcinn IVfiTTip 

1 V,Vj fcLlVJll i. 1 CLlllt* 


J J- mo ^ 




n ( Mr^ km 


N (cm -2 ) 




b A. Outflow 4a 


4.27 


.022 


.15 


1.4xl0 19 


11 


b B. Outflow 4b 


4.60 


.032 


.21 


1.5xl0 19 


13 


b C. Outflow In 


14.5 


.088 


.71 


4.8xl0 19 


66 


b D. Outflow 6/7 


4.45 


.045 


.30 


1.5xl0 19 


29 


r E. CO Region 3 


1.31 


.016 


.112 


4.3xl0 18 


8.5 


b F. Sh 2-233IR NE 


41.8 


.464 


3.72 


1.4xl0 20 


330 


m F. Sh 2-233IR NE 


132.9 


1.47 




4.4xl0 20 




1 F. Sh 2-233IR NE 


30.0 


.333 


2.03 


9.9xl0 19 


135 


b G. Outflowls 


14.6 


.064 


.48 


4.8xl0 19 


40 


r H. CO Region 2 


4.54 


.012 


.074 


1.5xl0 19 


5 


b I. Outflow 9 


6.33 


.039 


.39 


2.1 xlO 19 


43 


b J. CO Region 1 


3.61 


.015 


.12 


1.2xl0 19 


11 


r K. Red S 


5.26 


.051 


.34 


1.7xl0 19 


26 


b L. Blue S 


3.66 


.053 


.47 


1.2xl0 19 


47 


b 1' aperture 


15.1 


.96 


7.6 


5.0xl0 19 


670 


b 3' aperture 


2.7 


1.6 


12 


9.0xl0 18 


1000 


b 5' aperture 


1.7 


2.7 


20 


5.6xl0 18 


1600 


r 1' aperture 


11.8 


0.75 


4.7 


3.9xl0 19 


320 


r 3' aperture 


1.9 


1.1 


6.8 


6.2xl0 18 


460 


r 5' aperture 


0.96 


1.5 


10 


3.2xl0 18 


640 


b 1' 12 CO 2-1 


10.4 


.94 


7.1 


4.9 xlO 19 


590 


m l' 12 CO 2-1 


97.78 


8.83 




4.6 xlO 20 




r 1' 12 CO 2-1 


9.17 


0.83 


5.52 


4.3xl0 19 


430 


m l' 13 CO 2-1 


41.12 


211 




l.lxlO 22 




m l' C ls O 2-1 


5.31 


271 




1.4xl0 22 





a Unless labeled otherwise, regions are extracted from CO 3-2 data as shown in figure [6p 
b Blue integration over velocity range -34 to -21 km s -1 

c Apertures are centered on J(2000) = 05:39:11.238 +35:45:41.80 in Sh 2-233IR NE 
r Red integration over velocity range -13 to -4 km s -1 

m Middle range integration over -21 km s -1 to -13 km s - . Assumed not to be outflowing, so no momentum is 
computed 
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V LSR (km/s) 

Fig. 9.— CO spectra of inner 12" centered on Sh 2-233IR NE for all observed CO lines. The CO 
3-2 and 2-1 beams are not matched, but in both cases the area integrated over is 1-2 resolution 
elements across. The divisions demarcating the red and blue line wings are shown with vertical 
dashed lines at vlsr = — 21 and -12 km s _1 . 
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Fig. 10. — TripleSpec slits (blue) overlaid on the H2 image. The red boxes indicate the apertures 
extracted from those slits to fit and measure H2 properties. The apertures are also indicated in the 
position-velocity diagrams. 
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Fig. 11. — Position-velocity diagrams of the H2 2.1218 [im line in Outflows 1,2, 4, IR 6, and 
IR93/IR58. The velocity range is from -340 to 190 km s _1 . 

Table 3. TripleSpec fitted H2 outflow velocities 



Outflow Number Aperture Number a ^lsr (km s 1 ) b V£Sij(km s 1 ) 



1 


1 


-33.54 


(0 


15) 


-31.85 (0.32) 


1 


2 


-13.60 


(0 


57) 


-13.56 (0.96) 


1 


3 


-40.51 


(0 


41) 


-36.13 (0.81) 


1 


1 


-88.7 


(2 


8) 


-83.7 (7.9) 


2 


1 


-82.6 


(7 


6) 


-81 (21) 


2 


2 


-30.41 


(0 


57) 


-28.9 (1.4) 


2 


3 


-33.89 


(0 


62) 


-35.2 (3.7) 


4 


1 


-73.34 


(0 


48) 


-70.2 (1.1) 


4 


2 


-64.08 


(0 


61) 


-67.8 (2.2) 


IR6 


1 


-39.4 


(1 


6) 


-39.4 (4.2) 


IR93 


2 


-26.07 


(0 


43) 


-26.85 (0.97) 


IR93 


3 


-30.6 


(1 


5) 


-32.0 (2.5) 


IR93 


•1 


-29.14 


(0 


77) 


-30.3 (2.1) 


IR93 


6 


-47.7 


h 


9) 


-71 (37) 



a Measured from H 2 1-0 S(l) 2.1218 fim line 

b Measured from all detected H2 lines fit with model described in section [3.41 
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Table 4. Measured properties 


of H2 flows 






Outflow 


a Center 


b PA c Length 


d Source 


c Flow e Counterflow 


f Age 


s LOS 










Length Length 


(50 km s" 1 ) 


Velocity 


1 


05:39:13.023 +35:45:38.66 


-13.3 142.3" 


mm2? 


58 84.2 


1.4e4 




2 


05:39:13.058 +35:45:51.28 


-47.0 44.6" 


mmla 


44.6 


6.6e3 


Blue 


3 


05:39:12.48 +35:45:54.9 


-62 44" 


mm3? 


44 


6.5e3 


Red 


1 


ambiguous 


17.8-21.8 141-144" 


? 


141-144 


2.1c4 


Blue 


5 


05:39:12 +35:45:51 


170 38-48 


mm3? 


38-48 


6.5e3 


Blue 


6 


05:39:09.7 +35:45:17 


14.5 197 


Q5 


197 


2.9e4 


Blue 


8 


05:39:10.002 +35:45:10.87 


-154.6 105.5" 


IR41? 


54.7 52.9 


7.9e3 





a Midpoint of bipolar outflow if symmetric, position of jet source candidate if asymmetric 

b Position angle uncertainties are ~ 5° because they are not perfectly collimated, causing an ambiguity in their true directions. 
The exact angles used to draw vectors in figure [5] arc listed for reproducibility. 

c Total length of outflow on the sky, including counterflow 

d Candidate jet source object. Outflows 2 and 6 have clear associations, the others are weaker candidates. 

e Flow length is the distance from the CENTER position to the last H2 knot in the position angle direction as listed. 
Counterflow length is the distance from the CENTER position to the opposite far knot. 

f Timescale of jet assuming it is propagating at 50 km s — 1 , an effective lower limit to see H2 em ission. If two lengths are 
available, uses the longer of the two. These are lower limits to the true timescale (|Parker et al.lll99ll) . 

g The parity of the outflow along the line of sight. Outflow 1 and 8 have counterflows with parities as indicated in figure [2] 
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Table 5. Measured H2 line strengths 





1-0 S(0) 


1-0 S(l) 


1-0 S(2) 


1-0 S(3) 


1-0 S(6) 


1-0 S(7) 


1-0 S(8) 


1-0 S(9) 


1-0 Q(l) 


1-0 Q(2) 


1-0 Q(3) 


1-0 Q(4) 


aperture 


2.2233 


2.12183 


2.03376 


1.9S7S6 


1.7879S 


1.74803 


1.71466 


1.68772 


2.40659 


2.41344 


2.42373 


2.43749 


outflowlapl 


3.60E-15 


9.80E-15 


5.50E-15 


1.20E-14 


4.70E-15 


3.10E-1S 


8.60E-16 


1.10E-1S 


9.20E-15 


6.10E-1S 


1.10E-14 


6.90E-15 




( 2.4c-17) 


( 3.4c-17) 


( 6.8e-17) 


( 2c-16) 


( 2e-16) 


( 2.8c-17) 


( 2.8c-17) 


( 2.7e-17) 


( 1.4e-16) 


( 7.4c-17) 


( 8.8c-17) 


( 7.8c-17) 


outfiowlap2 


7.10E-16 


1.80E-15 


9.90E-16 


1.80E-15 










3.00E-15 


1.90E-15 


3.20E-15 


2.00E-15 




( 2.1c-17) 


( 2.7c-17) 


( 6.8e-17) 


( 1.7c-16) 










( 1.3e-16) 


( 4c-17) 


( 7.8c-17) 


( 3.8e-17) 


outnowlap3 


1.60E-15 


4.10E-15 


2.20E-15 


4.70E-15 




8.30E-16 






5.60E-15 


3.70E-15 


6.60E-15 


4.80E-15 




( 2.4c-17) 


( 3.4c-17) 


( 6.3e-17) 


( l.Sc-16) 




( 2.8c-17) 






( 1.4e-16) 


( 5.9e-17) 


( 8.2c-17) 


( 5.2e-17) 


outfiowlap4 




9.00E-16 


























( 3c-17) 






















outfiow2apl 




3.60E-16 


























( 1.5c-17) 






















outnow2ap2 


9.40E-16 


2.40E-15 


1.50E-15 


1.80E-15 




4.00E-16 






3.00E-15 




3.70E-15 






( 1.7c-17) 


( 2.2c-17) 


( B.8e-17) 


( l.le-16) 




( 2.3e-17) 






( 4.7c-17) 




( 7.9c-17) 




outfiow2ap3 


2.10E-1B 


1.90E-15 


1.80E-15 


2.20E-15 




6.70E-16 






5.70E-15 




7.30E-15 






( 1.7c-17) 


( 2.2c-17) 


( 5.1e-17) 


( 1.3c-16) 




( 2.9c-17) 






-8.00E-16 




-8.00E-16 




outfiow4apl 


5.50E-16 


2.00E-15 


8.50E-16 


2.00E-15 




9.40E-16 


1.90E-16 


3.40E-16 


1.40E-15 




1.40E-15 






( 2c-17) 


( 2c-17) 


( Be- 17) 


( 1.3e-16) 




( 2.8c-17) 


( 1.8c-17) 


( 2.3e-17) 


-4.00E-16 




( 6.9c-17) 




outfiow4ap2 


5.60E-16 


2.00E-15 


S.30E-16 


2.10E-1B 




5.80E-16 




1.10E-16 






2.00E-15 






( 2c-17) 


( 2.2c-17) 


( 2.4e-17) 


( 1.2c-16) 




( 2.3c-17) 




( 1.8c-17) 






-2.00E-16 




IR6apl 




1.10E-15 




9.30E-16 




4.30E-16 


















( 3e-17) 




( 1.4c-16) 




( 3.2e-17) 














IR93apl 




6.60E-15 




2.70E-15 














5.80E-15 








( 3.5e-17) 




( le-16) 














( 7.4c-17) 




IR93ap2 


4.40E-15 


6.60E-15 


3.90E-1S 


3.30E-15 




1.10E-1S 






7.60E-15 


S.10E-1S 


6.90E-15 


5.50E-15 




( 3.2c-17) 


( 3.7c-17) 


( 9.2e-17) 


( 1.4c-16) 




( 2.7c-17) 






( 8c- 17) 


( 5.2c-17) 


( 7.4c-17) 


( 6.1e-17) 


IR93ap3 


1.00E-15 


1.70E-15 




9.00E-16 










2.00E-15 


1.70E-15 


1.90E-15 






( 2.3e-17) 


( 3.6e-17) 




( 1.2e-16) 










( 8c- 17) 


( 3.8e-17) 


( 8.8c-17) 




IR93ap4 


2.60E-15 


3.70E-15 














4.30E-15 


3.50E-15 


4.70E-15 






( 3.2e-17) 


( 3.6e-17) 














( 8.5c-17) 


( 5.2c-17) 


( 7.4c-17) 




IR93ap5 




1.90E-15 




1.00E-15 






















(2.4e-17) 




(1.00e-16) 


















IR93ap6 




4.10E-16 


























( 3c-17) 























Note. — Fluxes arc in units erg s "'"cm 2 A 1 . Errors are listed on the second row for each aperture. Errors of (0) indicate that the line was 
detected, but that the fluxes should not be trusted because the background was probably oversubtractcd. 
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Table 5. Measured H2 line strengths (cont'd) 





2-1 S(l) 


2-1 S(3) 


3-2 S(3) 


3-2 S(4) 


4-3 S(5) 


[Fe II] 


[Fe II] 




2.24772 


2.07351 


2.2014 


2.12797 


2.20095 


1.6435 


1.2567 


outfiowlapl 


2.00B-15 


1.20E-15 


6.20E-16 


2.60E-16 


7.10E-16 


4.4c-15 


3.5e-15 




( 2.5c-17) 


( 3.5e-17) 


( 2.2o-17) 


( 1.6o-17) 


( 1.9c-17) 


( 7.9c-17) 


( 4c- 17) 


outfiowlap2 












6.7c-16 


3.1o-16 














( 7.8o-17) 


( 3.3c-17) 


outflowlap3 


9.90E-16 




5.70E-16 


2.50E-16 


6.40E-16 


1.3o-15 


5.7o-16 




( 2.6c-17) 




( 0) 


( 1.2o-17) 


( 0) 


( S.9c-17) 


( 4c- 17) 


outfiowlap4 
















outflow2apl 
















outflow2ap2 


6.40E-16 
















( 1.9c-17) 














outfiow2ap3 
















outfiow4apl 


4.30E-16 


4.20E-16 






1.60E-16 








( 1.9c-17) 


(0) 






(0) 






outflow4ap2 

















IR6apl 
IR93apl 

IR93ap2 3.S0E-15 - 3.10E-15 

( 2.2c-17) (0) 

IR93ap3 
IR93ap4 

IR93ap5 - 
IR93ap6 



Note. — Fluxes are in units erg s -1 cm _2 A -1 . Errors are listed on the second row for each aperture. Errors 
of (0) indicate that the line was detected, but that the fluxes should not be trusted because the background 
was probably ovcrsubtractcd. 
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3.5. Spectroscopic Results: Optical 

IR 6 and IR 41 (objects 1 and 6 in Figure [2]) both show Ha in emission. IR 41 is close to the 
reflection nebula in the southeast portion of IRAS 05358 and is probably the reflected star. The 
reflection nebula's spectrum is very similar to IR 41 's spectrum at Ha in both width and brightness 
(see Figure fl2l) . 




6545 6560 6575 
Wavelength (angstroms) 



Fig. 12. — A position-velocity diagram of IR 6 and 41 including the reflection nebula near IR 41 
(« 7.4M Q ). IR 6 shows a two-peaked Ha emission profile, but is the less massive (« 4.5M ) of the 
pair. The separation between the two sources is 55'/3, and each pixel is O'M. 

There are three components in the Ha profile of IR 6: a broad absorption feature seen far (~ 
400km s _1 from the line center) on the wings and two emission peaks. The peaks are separated by 
190 km s~ x and the blueshifted peak is weaker than the redshifted (Tabled]). The H/3 profile shows 
much deeper absorption and weaker emission but with similar characteristics. The presence of the 
Ha emission makes identification of the stellar type from the Ha line profile uncertain. The derived 



extin ction to IR 6 is at least Ay = 7 from an assumed Ha/H/3 ratio of 2.87 (jQsterbrock &: Ferland 



2006). The H/3 flux was measured from zero to the peaks of the emission profile and therefore 



probably overestimates the H/3 flux and underestimates the extinction. 



3.6. Radio Interferometry 



A point source was detected in the X, U, K and Q b and VLA maps with high significance at 
the same location as the X-band point source reported in lBeuther et al.l (|2007l ). Seven-parameter 
gaussians were fit to each image to measure the beam sizes and positions and flux densities. The 
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Table 6. IR 6 Deblended Profiles 





a Blue 


b Blue 


a Red 


b Red 


Absorption 


Gaussian / 


b Absorption 




Emission 


Wavelength 


Emission 


Wavelength 




Lorcntzian FWHM 


Wavelength 


Ha 


4.4X1CT 14 


6559.79 


1.3X10 -13 


6564.23 


-2.6xl0~ 14 


1.5 / 0.19 


6563.02 


H/3 


c 2.4X10- 14 


4857.68 


c 1.8xl0- 14 


4864.28 


d -4.6xl0~ 14 


0.17 / 16.5 


4861.91 



Note. — Measurements are made using a Voigt profile fit in the IRAF SPLOT task. 
a Flux measurements are in units of erg s _1 cm _2 A _1 

b Wavelengths are in Geocentric coordinates. Subtract 0.53A from Ho and 0.39A from H/3 to put in LSR coordinates. 
c H/3 emission was measured assuming a continuum of zero and therefore represents an upper limit in the H/3 emission 
d H/3 dcblcnding may contain systematic errors from a guessed subtraction of the H/3 emission 



Table 7. Lines observed in the optical spectra 



Source 



Ha 



H/3 [S II] 6716A [S II] 6731A [O I] 6300A [O I] 6363A [O I] 5577A 



Outflowl apl 
Outflowl ap2 
Ambient Medium - slit 1 
IR 41 nebula 
IR 41 
IR 6 



4.3xl0~ 16 

6561.49 
4.5xl0~ 16 

6561.22 
6.7xl0~ 17 5.3xl0" 18 

6562.87 4861.7 
2.6x10- 15 

6562.85 
6. 5x10" 15 

6562.9 

1.76X1Q- 13 a 4.1xl0~ 14 



5.7xl0~ 16 
6715.3 

4.5xl0~ 16 
6714.9 

l.OxlO" 17 
6716.7 



6.3xl0~ 16 
6729.6 

4.6xl0~ 16 
6729.3 

7.9xl0~ 18 
6731.2 



5.3xl0~ 16 

6299.7 
3.1xl0~ 16 

6299.4 
3.5xl0~ lfi 

6300.3 
4.4xl0- 16 

6300.3 
l.lxlO" 16 

6300.0 



1.8xl0~ 16 

6363.3 
1.2xl0- 16 

6363.2 
1.2xl0" 16 

6363.8 
1.9xl0- 16 

6363.9 
7xl0- 17 

6363.3 



4.8xl0~ 17 
5578.0 



Note. — Wavelengths listed are in A and are geocentric. To convert to LSR velocities, subtract 24.35 km s 1 . The ambient medium 
fluxes represent averages across the slit. Fluxes are in erg s -1 cm -2 A _1 . 

a H/3 measurement in IR 6 is an upper limit 
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measurements are listed in Table [8j The locations of the point source and the shape of the beams 
from the re-reduced X and Q band ima ges are d i splay ed in figure [T3j A Class II 6.7 GHz methanol 
maser was detected in IRAS 05358 bvlMentenl (|199ll ). It was observed with the European VLBI 
Network (EVN) by iMinier et al.l (|2000l ) and seen to consist of a linear string of maser spots that 
trace a probable disk in addition to maser spots scattered around a line perpendicular to the 
proposed disk (see Figure IT3"j) . The VLA source is more than a VL A beam away from the VLBI 
CH3OH maser disk identified by lMinier et al.l (|2000l ). It is to the southeast in the opposite direction 
of Outflow 2. Outflow 2 is at position angle -47°, while the disk is at PA 25°. The 8° difference from 
being perpendicular is well within the error associated with determining the angle of the outflow 
in this confused region, so the VLBI disk is a strong candidate for the source of Outflow 2. 



Diagram centered at mm1 a 5:39:1 3.08 +35:45:51 .3 




Delta RA (") 



Fig. 13. — A diagram of the region surrounding mmla from iBeuther et al.l (|2007l ). The ellipses 
are centered at the measured source cen ters and their sizes r epresent the beam sizes of the Plateau 
de Bure interferomete r at 1.2mm (blue. IBeuther et al.l 120071 ) . Gemini MICHELLE at 7.9/im (red, 
Longmore et al.ll2006l ). the VLA at 3.6cm (green), a nd the VLA a t 7mm (orange). The maser disk 
was measured with the European VLBI Network bv lMinier et al.l (|2000l ). so the size and direction 
of the disk are very well constrained. The black circle is centered on the pointing center of the VLBI 
observation and represents the absolute pointing uncertainty. The arrow pointing in the direction 
of Outflow 2 traces clumps along the outflow back to the mm emission region. The vector is not 
to scale - Outflow 2 is about 45" long. 



The astrometric uncertainty in VLA measurements are typically < 0.1". Different epochs 
of high-resolution X-band and Q-band data confirmed that the pointing accuracy is substantially 
better than 0.1" in this case. The V LBI absolute pointing uncertainty is reported to have an upper 
limit of ~ 0.03" (jMinier et al.ll2000l ). The separation between the VLA Q-band center and the VLBI 
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disk center is 0.22", whereas the separation between the combined X and Q band pointing centers 
is only 0.027", which can be viewed as a characteristic uncertainty. This is evidence for at least 
two distinct massive stars in a binary separated by ~ 400 AU. While the statistical significance of 
the binary separation is quite high using formal errors, the systematic errors cannot be constrained 
nearly as well. This object is a candidate binary system but is not yet confirmed. 



4. Analysis 

4.1. Near-Infrared Spectroscopic Extinction Measurements 

Extinction along a line of sight can be calculated using the 1-0 Q(3) / 1-0 S(l) line ratio. 



A\ = 1.09 



In 



Sy(S)/S u (Q) 

A ul (S)X Q /A ul (Q)X s 



A, 



-i 



(1) 



Because they are from the same upper state, their intensity ratio should be set by their Einstein A 
values times the relative energies of the transitions. However, as shown by iLuhman et al.l ([1998), 
narrow atmospheric absorption lines in the long wavelength portion of the K band, where the 
Q branch lines lie, can create a significant bias. Because the lines have not been corrected for 
atmospheric absorption, the Q branch fluxes should actually be lower limits. Since the 1-0 S(l) 



Table 8. VLA measurements near IRAS 05358 mmla 



Frequency 


Beam major / 


RA (error) 


Peak flux 


Map RMS 


Observed 


minor / PA 


Dec (error) 


(error) 


(mJy/beam) 


43.3 GHz 


0.022" / 0.029" / -10.4 ° 


05:39:13.065425 (0.000015) 


1.319 (0.027) 


0.179 






35:45:51.14732 (0.00031) 






22.5 GHz 


1.52" / 1.28" / 232 


05:39:13.05521 (0.0029) 


1.26 (0.04) 


0.091 






35:45:51.378 (0.046) 






15.0 GHz 


1.58" / 2.00" / ° 


05:39:13.062 (0.005) 


1.274 (0.065) 


0.124 






35:45:51.4 (0.1) 






8.4 GHz 


0.107" / 0.122" / 7.9 ° 


05:39:13.064548 (0.000036) 


0.506 (0.003) 


0.015 



35:45:51.170356 (0.000613) 



'Errors reported here are fit errors. Absolute flux calibration errors are negligible for the X-band 
data but are about equivalent to measurement errors for the K, and U bands and dominant in the 
Q band 
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transition at 2.1818 microns is affected very little by atmospheric absorption, and the exinction 
measured is proportional to the Q/S line ratio, the measured extinction should be a lower limit. 

The [Fe II] 1.6435 and 1.2567 \im lines were detected in Outflow 1, allowing for another 
direct measurement of the extinction. The measured ratio FR = 1.26/x m/L64jzm in Outflow 1 
was .8, while the true value is at least 1. 24 but may be as high as 1.49 (jSmith Hartiganl 120061 ; 
Luhman et al.l Il998l ; iGiannini et al.l 120081 ) . The extinction measured from this ratio ranges from 
Ay = 4.1 (FR = 1.24) to 5.8 (FR = 1.49). The S(l)/Q(3) ratio uncorrected for telluric absorption 
is .91, which yields an extinction lower limit of Ay = 18.7, is inconsistent with the measurement 
from [Fe II]. The Ha detection and H/3 upper limit give a lower limit on the extinction of Ay = 6.6, 
which is consistent with both of the other methods to within the calibration uncertainty. 

It is possible that the two measurements come from unresolved regions with different levels 
of extinction, though a factor of at least 3 change over an area ~ 100^4t7 far from the millimeter 
cores seems unlikely. A strong IR radiation field could plausibly change the line ratio from the 
expected Einstein A value. The question is not resolved but may be possible to address with 
near-IR observations of nearby bright HH flows with more careful atmospheric calibration. 



4.2. Optical Spectra 

4-2.1. Stellar Type 

IR 6 is suspected to be the source of the bright R~2 finger at PA ~ 15°. IR 6 is also a 24/im 
source and was detected by MSX (designation G173. 4956+02. 4218). We identify this star as a 
Herbig Ae/Be star. 



4-2.2. Density and Extinction Measurements 

The spectrum of knot Nl (the bow of Outflow 1) allowed a measurement of electron density 
in the shocks from the [S II] 6716/6731 line ratio . Densities were determined to be n = 700 cm" 3 
in the forward lump and n =500 cm~ 3 in the second lump. Ha, [N II] 6583, [O I] 6300, and 
[O III] 6363 were also detected, but no lines were detected in the blue portion of the spectrum 
presumably because of extinction. The measured velocities from [S II] are faster than the H2 
velocity measurements at about vlsr = — 68 db 5 km s^ 1 . 

There is also an ambient ionized medium that uniformly fills the slit with a [S II] -measured 
density n e = 120 cm -3 . Evidently, nearby massive stars are ionizing the low-density ISM located 
in front of IRAS 05358. This material is moving at velocity vlsr = — 7±5 km s _1 and is extincted 
by Ay = 1.5 as determined from Ha/H^ = 2.87 assuming the gas is at 10 4 K. 



-33- 



4.3. UCHII region measurement 



A uniform density, ideal HII region will have an intensity curve / 

-1.35 

K 



e TV ) where 



r 



1235 



HT 2 ( S 



GHz 



-2.1 



EM 



pc cur 



a(u, T) 



(2) 



following iRohlfs Wilson! (|2004l ) equation 9.35, where a{v, T) ~ 1 is a correction factor. By 
assuming an excitation temperature T ex = 8500K, blackbody with a turnover to an optically thin 
thermal source was fit to the centimeter SED. The turnover frequency from this fit is r =15.5 GHz, 
corresponding to an emission m easure EM = 7 . 4 x 10 8 pc cm -6 . This turnover frequency is lower 
than the ~ 35 GHz reported bv lBeuther et al.l ((2007). The turnover is clearly visible in the U, K 
and Q data points in figure PT4"1 



By assuming the X-band emission is optically thick, a source size can be derived. 



2r 



2k B T e 



-X 2 D 2 



1/2 



(3) 



where D is the distance to the source. Assuming a spherical UCHII region and a distance of 1.8 
kpc, the source has radius r =30 AU (for comparison, the Q band beam minor axis is ~90 AU, so 
the region could in principle be resolved by the VLA + Pie Town configuration). 

The measured density is n = (EM/r) 1 / 2 = 2.2 x 10 6 cm , with a correspo nding emitting 
mass M = n^m_H4/37rr 3 = 1.0 x 10 -6 M & using ^ = 1.4. Using Kurtz et al.1 (119941 ) equation 1, 



N Lyc = (8.04 x lO 46 ^ 1 )! 1 - - 85 



(4) 



the number of Lyman continuum photons per second required to maintain ioniza t ion is esti- 
mated to be N^yc = 5.9 x 10 44 , a factor of ~ 4 lower than me asured by 
closer to a B2 ZAMS star (~ 11M ) than Bl using Table 2 of |Panagia _ 



Beuther et al.l (|2007l ) and 
19731'! . If the star has not 



yet reached the main sequence, it could be significantly more massive (jHosokawa &: Omukaill2009l ) 
so our stellar mass estimate is a lower limit. 



The gravitational binding radius of a 11 M Q star is r g = 2GM/c 2 ~ 190 AU (the HII region 
is assumed to be supported entirely by thermal pressure, which provides an upper limit on the 
binding radius since turbulent pressure can exceed thermal pressure). The UCHII region radius 
of 30 AU is much smaller, indicating that, under the assumption of spherical symmetry, the HII 
region is bound. 



Leurini et al.l (120071 ) noted that the CH3CN line profile around this source could be fit with a 



binary system with separation < 1100AU and a total mass of 7-22 Mq. This is entirely consistent 
with our picture of a massive binary system with all Mq star in a UCHII region and another high 
mass star with a maser disk. 
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There are no other sources in the IRAS 05358 region to a 5a limit of 0.075 mJy in the X-band, 
which provides the strictest upper limit. From equation [31 this corresponds to an optically thick 
source size of 24 AU. The maser disk has a spatial extent of around 140 AU, so it is quite unlikely 
that either an undetected UCHII region or the observed UCHII are associated with the maser disk. 

Assuming the same turnover point for undetected sources, an upper limit is set on Ni yc for 
undetected sources: 

N Lyc = (8.04 x lO 46 *- 1 ) ( 2fc ^i,85 A ^^c) EM (5) 

Our 5a upper limit is Ni yc = 1.38 x 10 44 s _1 , indicating that any stars present must be a later 
class than B3, or lower than about 8 M@. For an emission measure as much as 3 times higher, the 
corresponding stellar mass would be less than 10 Mq. It is likely that no other massive stars have 
formed in Sh 2-233IR NE. 

After independent ly determining the b est-fit UCHII model to the VLA data, we included the 



PdBI data points from lBeuther et al.1 (120071 ) and fit a power-law to both data sets. If the emission 
measure was allowed to vary, the derived parameters were EM = 6.3 x 10 8 and j3 = 0.7. However, 
doing this visibly worsened the UCHII region fit without significantly improving the power-law 
fit, so the fit was repeated holding a fixed emission measure, yielding (3 = 0.8 (plotted in Figure 



114b). This power-law is much shallower than the (5 = 1.6 measured by iBeuther et al.l (120071 ) 
without access to the 44 GHz data point, and suggests that there is a significant population of 
large grains in source mmla. However, we caution that the fits were performed only accounting for 
statistical errors, not the significant and unknown systematic errors that are likely to be present 
in mm interferometric data. The PdBI beams are much larger than the VLA beams, so the larger 
beams could be systematically shifted up by including additional emission, which would reduce (5. 
Nonetheless, the new VLA data constrains the UCHII emission to contribute no more than 10% of 
the 3.1mm flux. 



4.4. Mass, Energy, and Momentum estimates from CO 

4-4- 1- Equations 
The column density for CO J=3-2 is estimated using the equation 



where A ui = A 32 = 2.5 x 10 _6 s _1 and A 2 \ = 6.9 x 10~ 7 s~ 1 ()Turner et al.lll977l ). the rotational 
constant B e = 57.64, 55.10, and 55.89 GHz for 12 CO, 13 CO, and C 18 respectively, r) mb = .68, and 
T ex is assumed to be 20K. The partition function is approximated as 



9 10 20 30 40 

Frequency (GHz) 



X 
13 



Frequency (GHz) 



Fig. 14. — (a) The HII region fit to measured X, K, U, and Q band data. Error bars represent 
statistical error in the flux measurement. The Q band error is dominated by flux calibration 
uncertainty (see Table [I]). The measured turnover is at 9.5 GHz, (b) A fit to both the VLA data 
presented in this paper and the (sub)mm points from iBeuther et al.l (|2007l ). The best fit spectral 
index for the dust emission is a = 2.8 (/? = 0.8), significantly lower than the a = 3.6 measured by 



Beuther et al.l (|2007l ) without the 0.7 mm data point. 
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which is valid when T ex >> hB e /kB ~ 2.8K. Equation [6] becomes 

iV H2 = (3.27 x I0 18 cm" 2 ) — / T* A {v)dv 
where the integrand is in units K km s _1 . The mass is then 

M = ii m H2 A iY Ha = 1.42 x 10" 5 ^^- J T* A (v)dv 



(8) 



(9) 



where A is the area in cm 2 , fi = 1.4 is a constant to account for the presence of helium, and again 
velocity is in km s _1 . 



{.2. CO J = 2-1 Isotopologue Comparison 



Thomas fe Fullerl ( 20081 ) observed C O in the J=2-l and 3-2 transitions each with a single 
pointing using the JCMT centered at J(2000) = 05:39:10.8 +35:45:16 and measured a column 
density iVjj 2 = 4.03 x 10 22 cm -2 . The peak column density is 1.7 x 10 22 cm~ 2 in 13 CO and 
2.2 x 10 22 in C 18 at J(2000) = 5:39:10.2 +35:45:26, which is reasonably consistent with the C 17 
measurement considering abundance uncertainties. The peaks of the integrated spectra for C 18 
and 13 CO are coincident, but the 12 CO integrated peak is at J(2000) = 5:39:12.6 +35:45:46 (Figure 
[3 discussed more in section 15 . 1|) . 

Measurements of the column density, mass, momentum, an d energy are performed as in Equa- 



tion[6l Assuming a CO/ CO ratio of 60 (ILucas fc Lisztlll998l ) and optically thin 3 CO, the mean 

The resulting total mas s of the central 



4.5 



10 21 cm- 



column density across the region is iVjj 

3' is about 320 M^, which i s subs tantially smaller than the 600 M measured bv lBeuther et al 



of 450 and 400 M from 



(|2002al ) and IZinchenko et al.l (119971). but it is nearly consistent with 870/um and NH3 estimates 



Mao Zend (|2004l ) and is within the systematic uncertainties of these 
measurements. Assuming C 18 is optically thin and the C 18 0/ 13 CO ratio is 10, the column den- 



sity is 5.2x10 cm and the mass is 360 M@, which is consistent with the CO measurements, 
indicating that optical depth effects are probably not responsible for the discrepancy with the dust 
mass estimate. 



4-4-3- CO Mass and Energy Measurements for Specific Outflows 

Table [2] lists measurements of mass and momentum in apertures shown in figure [6j Where red 
and blue masses are listed, there is an outflow in the red and blue along the line of sight. Where 
only one is listed, an excess to one side of the cloud rest velocity was detected and assumed to be 
accelerated gas from a protostellar outflow. Blue velocities are integrated from -33 to -21 km s _1 . 
Red velocities are integrated from -12 km s -1 to 1 km s _1 . All masses are computed assuming CO 
is optically thin in the outflow, which leads to a lower bound on the mass; 13 CO 2-1 was measured 
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to have an optical depth of 0.1 in 7 very high velocity outflow s in lChoi et al.l (|1993l ). so if a relative 
abundance 12 CO/ 13 CO= 60 is assumed (jLucas Lisztlll998l ). masses increase by a factor of 6. 



It is not possible to completely distinguish outflowing matter from the ambient medium. While 
the outflowing matter is generally at higher velocities, the outflow and ambient line profiles are 
blended. A uniform selection of high velocities was applied across the region but this may include 
some matter from the cloud, biasing the mass measurements upward. Outflows in the plane of the 
sky and low-velocity components of outflows will be blended with the cloud profile, which would 
lead to underestimates of the outflowing mass. The momentum measurements, however, should be 
more robust because they are weighted by velocity, and higher velocity material is more certainly 
outflowing. The momentum measurements are referenced to the central velocity of Sh 2-233IR NE, 
vlsr = -16.0 km s _1 . 



Discussion 



5.1. Outflow Mass and Momentum 



Beuther et all (|2002al ) reported a total outflowing mass of 20 M in Sh 2-233IR NE. We 
measure a significantly lower outflow mass of 2 Mq under the assumption that the gas is optically 
thin, but this assumption is not valid: a lower limit can be set from the weak 13 CO 2-1 outflow 
detection (lower l imit because not all of the outflowing material is detected) on the outflowing 
mass of ~ 4 Mq. IChoi et al.l ( 19931 ) measure an optical depth of 13 CO 2-1 r ~ 0.1 in 7 very high 
velocity outflows. Our 13 CO data sugges ts that the optical d epth is somewhat lower, r ks 0.07. 
The abundance 12 CO/ 13 CO= 60 is used (jLucas fc Lisztl Il998l ) to derive a total outflowing mass 
estimate M ~ 25 Mq. The total outflowing mass is therefore ~ 4% of the total cloud mass, though 
most of the outflowing material is coming from Sh 2-233IR NE, so as much as 13% of the material 
in Sh 2-233IR NE may be outflowing. 

The most prominent outflow in IRAS 05358, Outflow 1, is primarily along the plane of the sky, 
so the high velocity CO is likely associated with the other outflow s that have significant components 
along the line of sight. As pointed out in lBeuther et al.l (|2002al ). the integrated and peak CO are 
aligned with the main mm core. High-velocity H2 near the mm cores and the blueshifted outflows 
2 and 4 all suggest that there are many distinct outflows that together are responsible for the high 
velocity CO gas. 

The offset between the integrated 13 CO peak and 12 CO peak in the J=2-l integrated maps, 
which corresponds with an offset in the peak of the integrated CO 3-2 map and the peak temperature 
observed in CO 3-2, suggests that the gas mass is largely associated with Sh 2-233IR SW, but the 
outflowing gas is primarily associated with Sh 2-233IR NE. The integrated and maximum brightness 
temperatures in 13 CO and C 18 are also centered near Sh 2-233IR SW, which rules out optical depth 
as the cause of this offset. CO may be depleted in the dense mm cores, which would help account 
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for the lower mass estimate from CO isotopologues relative to dust mass and NH3. Alternately, the 
gas temperature in Sh 2-233IR SW may be significantly higher than in Sh 2-233IR NE except in 
the outflows, which are probably warm. In this case, the outflowing 12 CO enhances the integrated 
intensity because of its high temperature and reduced effective optical depth, but it does not set 
the peak brightness because of the low filling-factor of the high-temperature gas. 

Because the o utflows are seen in H2, which requires shock velocities ~ 30 km s _1 to be excited 
([Bally et al.ll2007l ). and because the association between the high- velocity CO and the plane-of-the- 
sky H2 is unclear, a velocity of 30 km s _1 is used when estimating the dynamical age. Assuming the 
outflow is about 0.5 pc long in one direction (e.g. Outflow 1), the dynamical age is 1.6xl0 4 years. 
Outflow 4, which is around 1 pc long, is also seen at a velocity of -70 km s _1 LSR, or about -50 
km s _1 with respect to the cloud, and therefore has a dynamic age 2x 10 4 years, which is consistent. 



5.2. Energy Injection / Ejection 

Using an assumed outflow lifetime of 5 x 10 3 years for v = 100 km s™ 1 as a lower limit (because 
the full extent of the flows is not necessarily observed) and 1 x 10 5 as an upper limit (for the CO 
velocities ~ 10 km s _1 and the longest ~ 1 pc flows), mechanical luminosities of the outflows 
L = E/t are derived. The summed mechanical luminosity of the outflows is compared to the 
turbulent decay luminosity within a 12", 1', and 5' radius centered on Sh 2-233IR NE in Table [9j 

The rate of turbulent decay can be estimated from the crossing-time of the region, L/v, where 
L is the length scale and v is the the typical turbulent velocity. On the largest (~ few pc) scales, the 
mechanical luminosity from high- velocity outflowing material is approximately capable of balancing 
turbulent decay and upholding the cloud against collapse. However, at the size scales of the Sh 2- 
233IR NE clump (~ 0.1 pc), turbulent decay occurs on more than an order of magnitude faster 
timescales than outflow energy injection. On the smallest scales, outflow energy can be lost from 
the cluster through collimated outflows, though wide-angle flows and wrapped up magnetic fields 



Table 9. Comparison of turbulent decay and outflow injection 



Radius (pc) 


hurt 3, (yr) 




Loutflows h (Lq) 


Binding Energy (erj 


;s) c Outflow Energy (ergs) 


Turbulent Energy (ergs) d 


0.10 


2xl0 4 


20 


0.03-0.6 


3.4 xlO 46 


3.5 XlO 44 


5.0X10 46 


0.52 


lxlO 5 


12 


0.6 - 9.4 


5.9 xlO 46 


5.9 xlO 45 


1.5X10 47 


2.62 


5xl0 5 


2.3 


1-22 


1.2 xlO 46 


1.4 xl0 4B 


1.5xl0 47 



a Masses are assumed to be 600 Mq for the 1' and 5' apertures, and 200Mq for the 12" aperture. 

b Outflow luminosities are given as a range with a lower limit L = E ou t/I0 5 jr and upper limit L = E ou t/5 X 10 3 yr, where E ou t is 
from Table [2] multiplied by 6 to correct for outflow opacity. 

c Binding energy is the order-of-magnitude estimate GM 2 /R 

d Turbulent energy is computed using the measured 5 km s — 1 line width as the turbulent velocity. 
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will not propagate outside of the core region. Once collimated flows impact the local interstellar 
medium in a bow shock, their energy and momentum are distributed more isotropically and again 
contribute to turbulence. The imbalance on a small size scale is consistent with the observed 
infall signature (Figure [9]) in the inner 12" around Sh 2-233IR NE and the lack of a similar profile 
elsewhere. 



5.3. Comparison to other clumps 



The classification scheme laid out in iKlein et al.l (|2005^ is used to ident ify Sh 2-233IR NE as 
a Protocluster and Sh 2-233IR SW as a Young Cluster. iMaurv et al.l (|2009l ) performed a similar 
analysis of the Early Protocluster NGC 2264-C. They also found that the outflow mechanical 
luminosity could provide the majority of t he turbulent ene r gy L± U rb ~ 1.2L Q within the protocluster 
in a radius of 0.7 pc with a mass 2300Mq. IWilliams et al.l (|2003l ) performed an outflow study of the 
OMC 2/3 region with radius 1.2 pc and mass 1100 M , which is also an Early Protocluster, and 
concluded that L tur & ~ Lfi ow ~ 1.3-Zvq. While all three regions have nearly the same turbulent decay 
luminosities and outflow mechanical luminosities, Sh 2-233IR NE in IRAS 05358 is significantly 
more compact and lower mass than the Early Clusters, and is the only one of the three that 
contains signatures of massive star formation. 



5.4. Surrounding Regions 

About 8' to the southeast of IRAS 05358 is another embedded star form ing region, G173.58+2.45 



Inter f erom etric and stellar populatio n studies have been performed by IShepherd &: Churchwel] 



(jl996l ) and lShepherd Watsonl ([2002]). The bipolar outflow detected in their interferometric maps 
is also cleanly resolved in our figure[6l In our wide-field H2 maps, there is a complex of outflows 
similar to that of IRAS 05358, but fainter. 

The large HII region Sharpless 231 to the northeast can be seen in the Ha image (figure [1]). 
The expanding HII region is pushing against the molecular ridge that includes IRAS 05358 and 
accelerating the CO gas in the blue direction (e.g. the northern blueshifted clumps in figures [6] 
and [5]). It can be seen from the IRAC 8fim data that UV radiation from the HII region reaches 
to the IRAS 05358 clusters. The expanding HII region's pressure on the molecular ridge may be 
responsible for triggering the collapse of IRAS 05358 and G173. The size gradient from S232 (~ 30' 
across) to S231 (~ 10') to S233 (~ 2 — 3') is suggestive of an age gradient assuming uniform HII 
region expansion velocities and a common distance. Investigation of this hypothesis will require 
detailed stellar population studies in the HII regions with proper regard for eliminating foreground 
and background sources. 
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5.5. Massive Star Binary 

Our identification of a probable massive star binary with an associated outflow con tributes to 



a very small sample of known maser disks with H2 emission perpendicular to the disk. |Pe Buizer 



(|2003l ) observed 28 methanol maser sources with linear distributions of maser spots in the H2 2.12 
\im line, and he identified only 2 sources with H2 emission perpendicular to the maser lines. None of 
the outflows identified in his survey were as collimated as Outflow 2, so the methanol disk / outflow 
combination presented here may be the most convincing association of a massive protostellar disk 
with a collimated outflow. 

The association of a massive star with an UCHII region and a methanol maser disk and the 
very small s i ze of the UCHII region both suggest that the massive stellar system is very young. 



Walsh et al.l (|1998l ) suggested that the development of a UCHII region leads to the destruction of 
maser emission regions. Their conclusion is consistent with our interpretation of mmla as a binary 
system. 



6. Summary &l Conclusion 

We have presented a multiwavelength study of the IRAS 05358 star forming region. IRAS 
05358 contains an embedded cluster of massive stars and is surrounded by outflows. The outflows 
were linked to probable sources and determined that at least one outflow is probably associated 
with a massive (~ 10M Q ) star. Added kinematic information and a wide field view of the infrared 
outflows has been used to develop a more complete picture of the region. 

• Sh 2-233IR NE is a Protocluster and Sh 2-233IR SW is a Young Cluster 

• Energy injection on the scales of IRAS 05358 can maintain turbulence, but on the small scales 
of the Sh 2-233IR NE protocluster, is inadequate by ~ 2 orders of magnitude. Sh 2-233IR NE 
is collapsing. 

• there are 11 candidate outflows, 7 of which have candidate counterflows, in the IRAS 05358 
complex 

• there is a probable massive binary with one member of mass 12 M & in mmla, and the other 
which is the source of Outflow 2 

• there are at least two moderate-mass (~5M Q ) young stars in IRAS 05358 

We have identified additional middle- and high-mass young stars with outflows, and presented 
a case for a high-mass binary system within the millimeter core mmla. 
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